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The t-J model in the spinless-fermion representation is studied. An effective Hamiltonian for 
the quasiparticles is derived using canonical transformation approach. It is shown that the rather 
simple form of the transformation generator allows to take into account effect of hole interaction 
with the short-range spin waves and to describe the single-hole groundstate. Obtained results are 
very close to ones of the self-consistent Born approximation. Further accounting for the long-range 
spin-wave interaction is possible on the perturbative basis. Both spin- wave exchange and an effective 
interaction due to minimization of the number of broken antiferromagnetic bonds are included in 
the effective quasiparticle interaction. Two-hole bound state problem is solved using Bethe-Salpeter 
equation. The only d-wave bound state is found to exist in the region of 1 < (t/J) < 5. Combined 
effect of the pairing interactions of both types is important to its formation. Discussion of the 
possible relation of the obtained results to the problem of superconductivity in real systems is 
presented. 
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I. INTRODUCTION 



The problem of the hole motion in an antiferromagnetic (AF) background of the local spins, originally arisen in 
connection with the study of the localized magnetic insulators, has received considerable attention since the 
discovery of the Cu02 based high-temperature superconductors. It is well established that at zero doping these 
materials are insulators with the long-range AF order, which is well described by the two-dimensional Heisenberg 
model. H The instability of long-range AF order under the small finite doping of carriers is due to the strong 
interaction of spins with mobile holes. The simplest model, which contains in itself this strong interaction is the 
t-J model. H Extensive studies of this model validity for the description of the real Cu02 plane result in the number 
of quantitative predictions for the range of parameters and in the set of possible t-J model generalizations. [7|-[TT| It 
is widely believed that the essential low-energy physics of the high-T c systems can be done using the pure t-J model 

in the standard notation of the constrained fermion creation (annihilation) operators c\ {ci_ a ) , (ij) denotes the nearest 
neighbor sites, Si is a local spin operator, iVj is the operator of the number of spins. For the consideration of the two- 
and many-hole problem one must include projection operators into J-term, which project out the unphysical states 
with both spin and hole at the same site. Such a procedure is described in Section V,B. Physically, i-term describes 
additional hole (singlet) hopping over the background of local spins, or, otherwise, the hopping of hole (vacancy) 
in the spin background. An important feature of this term is the absence of the double particle occupancy at site. 
Exclusion of doubly occupied states does not allow to apply mean field type approximations. 

The single-hole problem in the t-J model ([l]) has been extensively studied by the various analytical 
numerical p3|-p6[ techniques, which have provided the deep understanding of the character of the hole motion. For 
a review see, e.g., Refs. |2^£^|. Analytical results obtained within the self-consistent Born approximation (SCBA) 
2j| agree very well with the exact diagonalization studies on clusters, J28| variational, Jl6|-[l8| and other 
approaches. |19fl The main feature of hole motion revealed in these studies is the strong renormalization of the naive 
tight-binding result for the band due to hole "dressing" by the cloud of spin excitations, which leads to a narrow band 
(~ 2 J for t/J>l and ~ t 2 / J for t <C J) with minima at the ±(tt/2, ±7t/2) points on the boundary of the magnetic 
Brillouin zone (MBZ). 

The two-hole problem has received much attention due to the searching of possible pairing mechanisms. In spite of 
the large amount of work the full consensus on the problem of bound states in the t-J model is absent. There were a 
lot of works devoted to the study of the spin- fluctuation pairing and corresponding type of superconductivity [[30|-|34| . 
There is strong evidence that the long-range spin-wave exchange, which is the source of the dipolar interaction between 
holes |3^,|35|], can lead to the d- wave pairing in the t-J model. As it was established in Ref. |3(J the corresponding 
bound states are shallow and have the large size. Many efforts also have been aimed to the study of the t-J model 



1 



bound states originated from the fact that the two holes can minimize their energy by sharing the common link, that 
can lead to the picture of superconductivity by "preformed" pairs [ 571 . More specifically numerical works in exact 
diagonalization on small clusters and Monte-Carlo studies, which account for the latter interaction, provide negative 
energy of the bound state of the d x 2_ y 2 symmetry up to the values t/J ~ 3 — 5 [|38|-ffl||, which are relevant to the 
real compounds. Variational |^|, and some kind of quasiparticle calculations |Q yield the critical value of t ~ 2 J 
for this interaction, which is somewhat lower than the real one. Generally, there is no agreement on the energy of the 
groundstate of two holes and on their spatial correlation function |Q even between the similar approaches. 

In this paper we propose a canonical transformation approach to the t-J model problems, that allowed us to turn 
from the t-J model to an effective quasiparticle Hamiltonian, describing the "dressed" holes and their interaction of 
the "contact" type and via spin waves, and then to find the groundstate of two such quasiparticles. Both types of 
interactions of the most interest are accurately accounted for by our approach. In some sense, we use the ideas of the 
earlier works by Sushkov et al [pTp| j45| , |36t , where the same scheme was realized using quite different approach. 

To begin, let us describe the form of the Hamiltonian (jl]) we start with. The most popular analytical approach to the 
t-J model is the SCBA |^^3[20 2?J, which is based on the spinless-fermion representation for the fermion operators 
and Holstein-Primakoff |^|2(| or Dyson-Maleev 0] representation for the spin operators for the t-J model. Namely, 
this approach is applied to the spin-polaron Hamiltonian, which is followed from the t-J one (0) in the presence of 
the long-range AF order and in the linear spin-wave approximation 



where h^(h), aJ(a), are the spinless hole and magnon operators, respectively, 2Jcj q = 2J(1 — 7 q ) 1 ^ 2 is the spin-wave 
energy, Mk, q = 4(7k- q w q + 7k« q ), Uq,Vq are the Bogolubov canonical transformation parameters, 7k = (cosk x + 
cosfcj,)/2. The spinless-fermion representation fulfills the above mentioned constraint on double occupation exactly 
p0| and, therefore, the only approximation made is the spin- wave one. As it was recently shown in Ref. p9|, the two- 
loop corrections due to the higher order terms in (||) is analogous to the higher-order nonlinear spin-wave correction 
to the linear spin-wave theory, and have the same order of smallness. 

We will study this (g) version of the t-J model with the additional interaction terms from the projection operators 
in J-tcrm. In such a formulation (|^) the t-J problem is explicitly the problem with the extremely strong interaction. 
The problem of the interaction of fermion excitations with bosonic field and the resulting effective "dressing" of 
fermion by the virtual cloud of bosons is an old and well investigated problem, and the powerful approach to it is 
the canonical transformation one |47| . Therefore, one can hope that a canonical transformation will be found very 
helpful for the t-J model too. Briefly, we will show that the rather simple transformation, which takes into account 
the main effect of the strong interaction ~ t and allows to consider the rest of the interaction perturbatively, exists. 

To complete the consideration of the known facts about the Hamiltonian (J2j> let us note, that in the recent work 
by Reiter p8| an exact wave function of the single hole has been obtained within the SCBA 



hi + Mc, q Gk- q (£k - ^q)^L q a q + • ■ • + M k,qGk-q(£k - Wq ) X . . . (3) 

10} , 
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where Gk(w) is an exact hole Green function, is the hole energy. Since it is an exact wave function of the 
Hamiltonian (||) (within the SCBA) there is no interaction of the corresponding quasiparticle with spin-waves, i.e. the 
initially strong interaction is transformed exactly to the "dressing" of the bare hole, and to an effective interaction 
between such quasiparticles. Unfortunately, one cannot simply obtain the hole-hole interaction extracting from 
the wave function (^) as a Fermi operator and averaging H (Eq. (||)) over the two-hole wave function hUi^^O), since 
h, h\ defined by this way, do not obey the usual anticommutation relations. 

Briefly, we present an approximate solution of the diagonalization problem of the initial Hamiltonian (||). An 
effective Hamiltonian is formulated for the quasiparticles, which have the energy, bandwidth, and structure very close 
to SCBA ones. Further solving of the two-hole problem is straightforward. 

Described procedure is valid for the region < {t/J) < 5. We consider the region 1 < (tj J) < 5 as an actual one, 
since considering t/J model as a result of the simple Hubbard or many-band Hubbard model mapping, t/J parameter 
has the lower boundary t/J ~ 1 below that the mapping procedure is not valid. Moreover, t/J — 5 corresponds to 
U/t = 20, which is hardly realized in the real compounds. 

The paper is organized as follows. In Section II, we give the comparison of the lattice polaron problem with the 
spin-polaron one and write the general form of the transformed t-J Hamiltonian. In Sections III, IV, we apply the 
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proposed procedure to the Ising case, small t limit, and general case of the problem. Section V is devoted to the 
two-hole problem. 



II. CANONICAL TRANSFORMATION 



From the formal point of view the spin-polaron Hamiltonian (g) has the form that is very similar to one of the 
usual lattice polaron problem. We consider here the lattice polaron problem to compare these two models in detail, 
establish similarities and differences. 

The Frohlich Hamiltonian is 

H = E E ^ + E "q 5 ^ + E^ c k- q c k (fo q + 6_ q ), (4) 

k q k,q 

where c'(c) and b'(b) are the electron and phonon operators, E^ and f2 q are their energies, respectively. 7 q f2 q is the 
electron-phonon coupling. At the limit of the "static" electron (Eu — Eq) diagonalization of the Hamiltonian can be 
done exactly using the Lang-Firsov transformation Ejl : 

H eff = e- s H e s = H + [H, S] + i [[H , S] S] , 

S = -E^ c Uck(& q -6- q ). (5) 

k,q 

Only the first two commutators are not equal to zero in this limit. One easily obtains the effective Hamiltonian in 
the terms of the "dressed" electron 

Heff =U-E n M 2 ) E 4<* + E n « b l b « E ^ql7q| 2 4- q 4 +q Ck< C k. (6) 

\ q / k q k,k', q 

Thus the electron-phonon interaction term in Eq. (Q) results in the lowering of the electron energy (polaronic shift) 
and the direct nf n| interaction. For the mobile electron an infinite series of terms in Eq. (|J) arises. It is summed up 
and yields the effective hopping term describing the collective hopping process of the bare electron with the cloud of 
phonons. It was shown that this "dressing" leads to the exponentially narrow effective band jl7|]. The emitting of the 
phonons (multiple phonon processes) can be considered as the perturbation. In the basis of this approach there is a 
clear physical idea that the presence of the electron at the lattice site leads to the change of the equilibrium positions 
of the surrounding ions, and the new eigenfunction of phonons is a coherent state. 

There are two main differences between the phonon and magnetic polaron problems. The first one is the principal 
absence of the "bare" dispersion in the Hamiltonian (||), i.e. its hopping term is more likely the vertex than the usual 
tight-binding hopping integral | p7[ . The second one is the nonlocal character of the hole-spin interaction, i.e. emission 
(absorption) of magnon can be done only by hopping. Because of that there is no "static" limit of the problem even 
if t -C J, and the evident a priori ideas about the structure of spin cloud around the hole are absent. 

Nevertheless, the existing knowledge of the character of hole motion in AF background can help one to succeed 
in turning to an effective model, which is much more appropriate to study than the initial one. Firstly, at the Ising 
background the groundstate of the hole is the localized magnetic polaron, which is formed by self-retraceable motion 
of the hole. At the Neel background there is the similar situation, i.e. spin waves in the virtual spin cloud around the 
hole are absorbed exactly in the reversed order than they were emitted. The contribution of the processes beyond 
these retraceable paths (or SCBA) approximation is found to be extraordinary small. Secondly, it was found in the 
number of works that taking into account the hole "dressing" even by the single spin wave already provides results, 



which are close to the exact ones |17|, 18 . Namely, the bottom of the band, hole minima locations, and width of the 
band are determined with a sufficient accuracy [ fi"7j ]. Therefore, it shows that the main contribution to the polaron 
well formation for the actual range of (t/ J) < 5 is made by the "one-string" component of the hole wave- function (p)) . 
Authors of some SCBA works also successfully used this approximation for some other t-J model studies p8[ , p9| . 
These are the reasons to hope that relatively simple transformation in the spirit of Lang-Firsov one can be used to 
get the effective model, which accounts for the main polaron effect (of the order of t) in the hole energy and hole-hole 
interaction, whereas the other included terms allow to apply the perturbation theory. 
We propose the general form of the generator of such a transformation: 

s=E^(4>4-h4 (7) 

k,q 
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where /Xk. q is the parameter of the transformation. It is natural to require for /Xk, q to obey the same symmetry 
properties as the kinematic factor Mk iq of the i-term in the Hamiltonian (0). Note that Mk iq is odd with respect to 
the transformations Mk, q = — Mk+Q, q = — A/k, q +Q, here Q = (77, it). So, without loss of generality one can rewrite 
fJ-k.q = ik, q -Mk, q , where /k, q is even under mentioned symmetry transformations. 
Main terms of the spin-hole Hamiltonian (0), which we wish to decouple are: 

H = tJ2 M k,q + H.C.) + 2 J ^ CJqflqflq, (8) 

k,q q 

after the applying of the transformation (^) they provide transformed Hamiltonian 

H tr = E E k hlh k + 2jJ2 w q a q a q +i E ^^ q 4- q 4'+ q ^k + tj^ ^k, q A/k, q (4-^4 + H.c.) 

k q k,k', q k, q 

+t ^(k, q, q') (4_ q+q ,ft k a q a q , + H.c.) + i E ^"-(k, q, q')4- q+q ^ka q a q , + • • • , (9) 

k, q , q ' k, q , q ' 

where we omit terms with six and more fermion operators. General expressions for the hole energy E k , hole-magnon 
formfactor Fk. q (up to the sixth order of the transformation), hole-hole vertex V^, (up to the fourth order) and 
other vertices are given in Appendix A. There is the freedom in the choosing of the parameter of the transformation 
/k. q - By analogy with an exact single-hole wave function (||) one can propose the following form of the parameter 
/k, q = /kG(k — q, Ek — 2Jojq), where G(k, u>) is the Green function of the dressed hole. As a result one obtains the 
very complex self-consistent integral equation on E k and self-energy Ek(cc), which is hardly soluble. There are some 
other ways of choosing the transformation parameter (TP), which allow to avoid self-consistency in equations and 
technically are more advantageous. The simplest one is to neglect q-dependence in /k. q => /k, and then to determine 
it as a solution of an equation obtained from the requirement of the minimum of the hole energy, or, for example, 
equality to zero of the hole-magnon formfactor. We will discuss this approach in the next two Sections. Here we claim 
that for the rather general form of TP one can restrict oneself by the first four terms in the transformed Hamiltonian 
@, namely 

Heff = E E k hih k + 2jj2 ^q« q « q + 1 E <%,A-X'+A h * + 1 E F ^ q M k , q (&L q fck4 + H - c -) ( 10 ) 

k q k,k', q k, q 

keeping in mind that the transformed "additional" interactions (J-term) is included in V£^, . Moreover, resulting 
effective hole-magnon vertex is perturbative, i.e. its second-order addition to the energy 5E k is small. Importance of 
this vertex for the two-hole problem we will discuss in the Section V. 



III. ISING LIMIT 

Let us start the general consideration of our approach from the Ising case. As it was noted in Ref. p0|| treating the 
t-J model in the Ising limit within the linear spin- wave approximation remains the physics of the problem essentially 
unchanged. Moreover, it was shown |50|| that the spin-wave formalism provides exactly the same result as one of the 
SCBA. 

One can get the Ising limit of the t and J terms from the general spin-hole Hamiltonian Eq. (||) using the 
momentum- independence of w q at the Ising background and the absence of the spin fluctuations (u q = 1 and v q = 0): 

H = t e ( h l-^ h A + H - c -) + 2 J E a A ( n ) 

k,q q 
With = 47k-q, 

the additional terms of the interaction Hamiltonian (from J-term) can be considered independently and we will include 
them later. 

Following the analogy with the Lang-Firsov transformation we turn to the effective Hamiltonian with the help of 
the transformation: 

H ef f = e- s He s = H + [H, S] + i [[H, S]S] + . . . , 

S = f E M k, q (ftl- q ftk< ~ H.c.) , (12) 

k, q 
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where the generator of the transformation simply reproduces the kinematic structure of the hopping Hamiltonian and 
involves the single free parameter /. It is natural for parameter of the transformation / to be k-independent, since 
the energy of the hole for the Ising problem does not depend on k. Using the evident relation [Hj,S] — f(2J/t)H t 
one can make the effective Hamiltonian calculation and get: 



Heff =E h J2 hlhk + 2JJ2 «q«q + t E ^k', q 4- q 4' +q ^/ik 
k q k,k',q 

+i^X>k, q (4- q >*k4 + H - c -) + t E ^ q V l k-q+q'^ q a q ' + . . . 
k, q k, q , q ' 

with the one-hole energy, hole-magnon formfactor, hole-hole vertex, and hole-two magnon vertex: 

2J 



(13) 



F 



k,k',q 



8t (^/-|/ 3 + 

i-4/ 2 + ^(/-|/ 3 ) 



(14) 



q,q 



[Ml .Ml 



k-q.-q 



Ml,,„ Ml, 



i • jf If (i - ^/ ; 



8 f 3 I i 

1 k, q 1 "k— q ,- q T~ Ju k'+ q , q JM k',- q j 3J I 



J 

2T 



Vfc/ - (Mi q Mi 



q"k q <| q' "k q',— ij'M^q/q) ■ 2/ [ 1 + / 



,/ „i - ML 



up to the fourth order of transformation for the hole energy and formfactor of the spin-hole vertex and for the hole- hole 
interaction and higher vertices. The first peculiar feature of the Ising case appeares here. Minimization of the energy 
provides an equation on /: 



^~l-4/ 2 + (2JA) (/-§/ 3 )=0 



(15) 



that coincides with equality to zero of the hole-magnon formfactor F. This is closely connected to the fact that the 
each act of emission or absorption of magnon is due to the hole hopping, and the polaron is created by moving on 
the self-retraceable paths. The role of the so called Trugman processes [^5) among the other fifth and sixth order 
contributions was found small (~ (1/8) of relative magnitude) and they remain all results essentially unchanged. 
The next simplifying fact is the absence of the two-magnon vertices with the h^ha^a^ (aa) terms. It means that at 
zero temperature there are no contributions of the hole-two-magnon part into the self-energy and to the hole-hole 
vertex, and hence, the h^ha^a can be omitted in the effective model. Thus, after energy minimization the effective 
Hamiltonian has the form, which is very similar to the lattice polaron one: 



H, 



eff 



= E h Y, h l h k + 2 JE fl q fl q + 1 E 



V hh 



(16) 



k,k', q 



here the ener gy a nd hole-hole vertex are given by Eq. (|14j) with / obtained from Eq. (|15| 
Eqs. (|l4|), ( |l5|) for the hole energy show that: 



f 2J' 



t/J ^ 



(17) 



demonstrating the perturbative nature of our approach at small t/J and some kind of l/z expansion at large t/J. 

An exact result for the single-hole energy on the Ising background was obtained in Ref. |5(J in the form of the 
difference equation. Also, there is an analytical solution of this equation in the f/J > 1 limit first proposed by 
Bulaevskii, Nagaev, and Khomskii H: 



E = -2yfzt -2J + 2. 34(2 J) 2/3 (V^) 



a/3 



(18) 
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Figure [I] presents numerical solution of the exact equation j5(| (bold solid curve) and approximate solution ( Jlq ) 
(dashed curve) together with our results Eq. (|l4|). Upper and lower curves correspond to the hole energy calculated 
up to the fourth and sixth order of the transformation, respectively. Figure || shows the components of magnetic 
polaron wave function for the exact solution (bold curves) and ones of our canonically transformed polaron: 

kt\0) = e - s rf\0) (19) 

These figures demonstrate that our single-hole energy is very close to the exact one, the higher orders for the actual 
tj J region play the role of corrections, and the components of exact and our polaron wave functions are close to each 
other. Therefore, one can hope that the consideration of the interaction of our quasiparticles will reveal the properties 
of exact eigenstates of the t-J Hamiltonian. 

IV. NEEL CASE. 
A. Small t limit. 

Let us consider first the small tj J limit of the model. Perturbation theory over this parameter works very well 
and the most results can be obtained analytically. Our transformation procedure in this limit also has a perturbative 
sense and it is possible to compare our results with ones of the usual perturbation theory. It is also useful to consider 
this limit for the demonstration of some details of our approach. 

As it was mentioned in Section II, there are some evident forms of the TP /k, q - The first one we wish to study is 
/ kq =>■ fa. The second one is / kq =>■ fa/^Oq, which can be considered as the simplification of one proposed in Section 
II /k,q =>■ /kG(k — q, — 2 Jujq) in the small t limit when <§; 2 Jw q . Technical advantage of both cases for any t 
is that k and q dependent parts are separable and all arisen integrals can be reduced to the several functions. 

It is interesting to clarify the physical meaning of both transformations. Let us show the polaron wave functions: 

1st case ^kJO) = 
2nd case ^ k |0) = 

Since at q — > -M^q — > the admixture of the long-range magnons in the wave function of the "first case" polaron 
is small, i.e. this transformation corresponds to taking into account short-range spin- wave "dressing" of the hole. 
Comparing the Reiter's wave function (||) with the "second-case" one suggests their identity for this limit. 
Transformed Hamiltonian (P) for the small t limit takes the form: 

q 

+f^M k!q (l + y/ MWq ) (4_ q fc k a q + H.c.) (21) 

k,q ^ ' 

+ t M k , q M k / +q , q / k , +q , q (l + y/k.q^q) ('lk-q4'+<A<^k + H.C.) + O(^j) , 

k,k',q ^ ' 

here we omit hole-two-magnon vertices (0(t ■ / J )) since they contribute in the energy and scattering amplitude only 
in the second order (0(t 4 / J 4 )). The hole energy for both types of the transformation: 

{J ~ 2 
Ik fa + —fafa 1st case 
( J \ ( 22 ) 
/ k 7 k (l + -/ k J 2nd case, 

with useful notations fa = ^ q M^ q , fa = X^q^kq^q* ^ k = E q ^k q / W q- Minimization of the energy averaging 
over the band provides the following parameters of the transformations and formfactors: 



;/kE M ^ ^ + /kE M ^U<+o((^) 3 ) io> 



K + fa 



Ah 



q 



L0 o 



q°q 



0((t/jf 



|0) . 



(20) 



H, 



eff 



E 



k.q 



1 + 4 



k,qWq 



G 




F k , q =<! * "4 1StCaSe , (23) 
2nd case 



and the effective Hamiltonians 
a n 



H eff = ~ 



1 fc k + 2 J £ ^ q a q a q + * ^ M k , q (l - (^_ q /» k a q + H.c 

k "'k q k,q ^/ 



57 M k,q^k'+q,q~ I * - ^q~~- I 

AJ k,k',q 7 k'+q V 27 k / 



(/>L q 4' +q V/ik + H.c.) (24) 



|rE^k + 2j£^ q a ' E Mk,q ^ k ' +q ' q (C q ^ +q ^fek + H.c. 

k q k,k',q q 

The hole energies have the well known shape with the minima at ±(7r/2, ±7r/2) points and large effective mass along 
MBZ boundary. Note, that the "second-case" transformation, which is the Lang-Firsov one by the construction, 
diagonalizes starting Hamiltonian exactly. Evidently, the "second-case" results for the hole energy and hole-hole vertex 
coincide with these values obtained by the usual perturbation theory. The first case correspond to the separation 
of the scales and implies the perturbative account for the interaction with the long-range spin waves. Clearly, that 
accounting for the second-order contribution of the remaindering part of the hole-magnon vertex (iJ e 1 // @) into the 
hole energy and hole-hole vertex reproduces results of the simple perturbation theory. 



B. General case. 



At the first glance the use of the second-case approach at arbitrary t is more appropriate, since one can hope for 
better account of the long-range spin waves into the direct hole-hole interaction and for small effective formfactor, 



which will allow to omit hole-magnon vertex at all and to get the simplest (24) effective model. At the large t/J (> 1) 
it is not correct, i.e. this transformation cannot reduce the initial hole-magnon vertex exactly to zero and about a half 
of the long-range spin-wave interaction remains untransformed. Therefore, the rest of the hole-magnon interaction 
cannot be neglected and the general form of the second-case effective Hamiltonian coincides with the first-case one. 
This is due to the difference of the Green function from the simple (— l/w q ) form at t/J > 1. 

At the same time, the "first-case" approach is better in a certain sense, since the resulting long-range part of the 
hole-hole interaction is separated from the short-range one, i.e. q at q — » tends to zero and the whole long-range 
part of interaction is from the rest of the single-magnon exchange. We focus on the long-range part of interaction 
since, as it was found earlier [Q, it is the key pairing interaction for the d x 2_ y 2 two- hole bound state. Thus, for the 
consideration of the two-hole problem correct account of the long-range magnon exchange is very important. In the 
situation when we technically can try the energy and formfactor up to the sixth order of the transformation, while 
the hole-hole vertex up to the fourth, the "first case" approach (/k, q => /k) becomes preferable. 

Shortly, we have performed the transformation using both approaches without significant differences in results, 
except the long-range magnon exchange amplitude, which 1.5 times less for the second case. Further we will discuss 
only the "first-case" results. 

Thus, we transform the initial Hamiltonian H = H t -j to an effective one H e ff 

H eff = e- s He s = H+[H,S] + j[[H, S]S] + ..., 

with (25) 
S = J2 /kA4, q (4-<A< - H.c.) , 

k,q 

which has the form of one in Eq. (|^) with parameters presented in Appendix A. 

For the sake of simplifying the consideration of the technical details of our approach let us omit for a moment the 
fifth and sixth order terms in the general formulae for the hole energy (see Appendix A). Minimization of the average 
energy by variation over the TP /k leads to the following integral equation: 
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r*k + g/k^k (26) 



= 



- f/k-fk-fk - ^/k-Ffc + |/k^kik ^ 3^/k-Wk 

where we use auxiliary functions i 7 ^ and Y"k, ^k given in Appendix A and notations introduced in the previous Section: 
I k = J2 M k,q = A o + ^i,o7k + ^i,i(7 k ) 2 , h. = M k, q ^q = 4> + A,o7k + A,i(7k f > ( 27 ) 



q 



with the short-hand notations 7k = (cos(k x ) + cos(fc 3/ ))/2, 7 k = (cos(fe x ) — cos(fc y ))/2 and numbers: 

A Q = 4.3273 , Ai j0 = -3.5907 , A^i = -0.2453 , (28) 
I = 3.7828 , A li0 = -3.2808 , A 1>t = -0.1496 . 

We use the following way of solving such a type of integral equation (^6|) . According to the discussion of the symmetry 



properties of the TP /k given in Section II and using the form of Eq. (J26J) it is evident that /k = /-k = /k+(7r,7r) = 
f(k x «-> k y ), and hence, it can be expressed as a power series in cos^) 2 , cos(fcj,) 2 , and cos(k x ) cos(k y ) , or more 
convenient 

oo 

/k= £ C„, m 7k ( "~ m) (7k ) 2m = Co,o + C lfil l + C M ( 7k ) 2 + ■■■ (29) 

then, substituting this form of /k in expressions for the auxiliary functions F£ (Appendix A), one gets an infinite 
number of integrals of the type ^ q q7k^q~"^ (7k-q) 2 ™) > eacn °f them is a finite series in 7^ (7k_ q ) 2 of the power 
(n + 2). Cutting /k and all other series at the finite power n one gets from Eq. (26) a set of (n + l)(n + 2)/2 nonlinear 
algebraic equations on Ci.j (i < n). Thus, the integral equation (|2^) is transformed to the set of algebraic equations, 
which is much easier to solve. Keeping in mind 1/z character of the energy expansion one can hope that only a few 
first terms are important, and the role of the higher orders is insignificant. 

We solved these systems of equations numerically for the particular values of0<t/J<5, and found that extension 
of the series in Eqs. (p6|), ( p9| ) from n — 3 (cos 6 , 10 equations) to n — 5 (cos 10 , 21 equations) changes results for 
the parameter /k, energy and formfactor (coefficients in their series) for the relative value less than 0.5%. Note, 
that including of the fifth and sixth order terms into expression of the energy (see Appendix A) change results for 
approximately 10%. In all further calculations we used the largest (n — 5) set of equations. 

With the solution for /k °f such a high accuracy in hand one can get explicit expressions for the energy, formfactor, 
hole-hole, and hole-two-magnon vertices for the effective Hamiltonian Eq. ([)]). Evidently, as in the small t limit, the 
hole energy has the shape with the minima at ±(7r/2, ±7r/2) points and large effective mass along MBZ boundary, 
and obeys the symmetry property — E^^^y 

The next step of our consideration is to prove the negligible role of the hole-two-magnon vertices and the perturbative 
character of the renormalized hole-magnon one. We have calculated the second-order corrections to the single-hole 
energy from the one-magnon and two-magnon self-energy diagrams for the various t/J. Briefly, correction to the 
depth of the band from the rest of the hole-magnon vertex is no more than 10%, while correction from the hole-two- 
magnon vertex (|^) is of the next order of smallness. Namely, for t/J = 3 E^/ 2 ,tt/2) = ~ 2.22t, 8E^ — — 0.15£, and 
5E^ = -0.02i. For the correction to the effective hole-hole vertex the relative contribution of the hole-two-magnon 
exchange is even smaller. Single magnon exchange is really negligible for the large transfer momentum, but is very 
important for the small one, indeed it has a "quasi-singular" form ~ t(q x + q y ) 2 /q 2 for the holes near the bottom 
of the band. As it was discussed above, this separation of the contributions to the effective hole-hole interaction in 
the momentum space leads from the construction of the transformation. Note also, that the two-magnon exchange 
cannot provide the singular interaction anywhere. 

Therefore, our analysis has shown the negligible role of the higher magnon vertices and proved that the using of the 
rather general type of the transformation leads to the transfer of the initially strong hole-magnon interaction mainly 
into the single-hole dispersion and hole- hole interaction. Thus, for a wide region of t/J and with the high level of 
accuracy one can restrict oneself by consideration of the effective Hamiltonian ( fL0| ) 

^e//=£Sk4^ + 2j£a; q a q a q + ^^ k , q Mk )q (4_ q /i k a q + H.c.) +t £ ^ q 4- q 4'+c^k'/>k , (30) 

k q k,q k.k'.q 
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with all quantities denned in the Appendix A, expressed through /ik.q = /k-^k.q, where /k is defined from the integral 
equation of the type of Eq. (^) . 

Figures |^-|| represent our results for the bottom and width of the single-hole band, and for the weights of the 
components of the magnetic polaron wave function together with ones from the SCBA calculations from Refs. [B0[ , 
p9| . From these figures we notice that the bottoms slightly differ for the all region of t/J, whereas the bandwidth 
difference is absent at small t and tends to rise at larger (t/J) > 5. The gap between the curve corresponding to our 
result and SCBA one in Fig. || is obviously due to the absence of the long-range magnon contribution in the former. 
Taking into account the second-order correction to the energy from the long-range magnon virtual emission-absorption 



provides the depth of the band Et- 



and difference E, 



(tt.O) 



E, 



(7T/2,7T/2) 



which almost coincide with the SCBA 



ones. Surprisingly, at small t/J (< 1, > 0.5) results for bandwidth are very close to each other (Fig. 

Figure || shows t/J dependence of the hole-magnon formfactor q at different momenta k and q. Following the 
general procedure of expansion in a series one can get the general form of the formfactor 



-Pk,q — 



E 



[(w n 



\ 2(r 
n,m,z,j/ 7k 



(7k) 2m 7^ q j) (7 k _ q ) 



2./ 



(31) 



= W ,o,o,o + Wi, ,o,o7k + W ,o,i,o7k-q + • • • + w q (W£o,o,o + VF 1 w 0j0:0 7k + W " 04i0 7k_ q + . . .) 



Thus, our formfactor has strong k and q-dependence at any t/J that differs from the result of Ref. ^] where the 
some kind of variational approach has been used. Namely, at t/J = 3 Fk,q changes from 0.4 at small q to 0.0 at large 
(~ 7r/2,7r/2) one. 



V. TWO-HOLE PROBLEM. 
A. Two-sublattice representation. 

Because of the AF long-range order there are two types of fermion and boson excitations associated with two 
sublattices. For considering the one-particle subspace it is of no importance whether one has the model with two 
different sorts of quasiparticles having the same properties (energy, wave function etc.), or the model with one type 
of quasiparticles. In the above consideration we used the latter for the sake of simplifying the notations. One can 
easily prove the formal equivalence of these approaches. For the two-sublattice representation there are two types of 
holes and magnons both defined in the first magnetic Brillouin zone, whereas for the one-sublattice representations 
holes and magnons are defined inside the full Brillouin zone, because of that E^ — E^ + ^ ^y 

For the calculation of the correlation function ]49| , consideration of the hole- hole interaction |51{] , or some other 
calculations into two-hole subspace one should turn back to the two-sublattice representation. It is convenient to do 
it using the following expressions of the and a q operators through the linear combinations of the new operators 

h k = (/ k + 9k)/V2 , h k+(w>w) = (/ k - 9k )/\/2 , (32) 
Oq = («q + /3q)/ \/2 , dq+^.Tr) = (aq - /3q)/ v2 , 

where /k and correspond to the fermionic excitations at the A and B sublattice, respectively. a q and /3 q are 
the two type of Bogolubov spin-wave excitations. Transition to these new variables for the hole-magnon part of the 
effective Hamiltonian is straightforward 

H eff => t E ^,qM k ,q (/ k _ q 5k/3q + ^-q/k^ + H.C.) , (33) 
k,q 

where summation is performed over the MBZ. 

Now we should obtain h,h—* f,g transformation of the hole- hole interaction part of the effective Hamiltonian (pSC]). 
Rewriting the hh interaction term ( |30|) in the form where summation is produced over MBZ leads to 

H# f => H*° + Hff +H^=tJ2 <k',q/Lq9k' + q^/k + * £ (^,q/Lq/^ + q/k' h + (/ "> <?)) • (34) 

k,k',q k,k',q 

Thus, there are three different parts in the H^ff, which correspond to the interaction between holes at the different 
sublattices (Jg-part) and at the same one (//- and gg-parts). The first terms for the interaction of the particles at 







the same sublattice arise in the third order of the transformation and physically correspond to the process shown in 
the diagram Fig. Qa. At the Ising background it is equivalent to the four-hopping process of the two holes staying 
at the apart corners of the minimal lattice square, when they move in the clockwise or anticlockwise direction doing 
the half of a turn. Generally, // or gg interaction has no some important features of the fg one. Namely, there are 
no singularities in their long-range interaction as well as the gain of the energy due to reducing of the number of 
broken AF bonds is absent. These physical reasons were checked earlier (i^] and it was found that the bound states 
are absent for the particles at the same sublattice in the region of (t/J) > 1. Thus, we will restrict ourselves by 
consideration of the interaction of the particles at the different sublattices. 

For obtaining fg interaction from hh one accurate consideration of the parity of the vertex V under the trans- 
formation R = k(k') — > k(k') + (7r, 7r) is required. There are two contributions of different parity (R = =f) in the 
effective fg interaction. Their diagrammatic analogues are presented in Figs. 0b and ^c, respectively. The first one 
is due to the one-magnon exchange and by its origin it is of the "exchange" type (V^, ). The second one is due to the 
two-magnon exchange and the contact interaction (additions to J-term) and it is of "direct" type (V^). Obviously, 
these contributions enter in the fg vertex with different sign 

<k', q = (-K'^(k, k', q) V^(k', k, -q) + V^(k, k', q) + V^(k', k, -q)) . (35) 

Note here, that the first non zero correction beyond the ladder approximation for the hole-hole (fg) scattering arises 
only in the sixth order over t (see Fig. ^i), that is the same order as one of the Trugman processes for the single 
hole movement and, as it seems, it has the similar reasons in the basis. Therefore, keeping in mind the negligible 
role of these non SCBA contributions to the hole energy, one can hope that the role of the diagram in Fig. |?]d in the 
hole-hole interaction can be omitted and the ladder approximation will work well even for the initial (untransformed) 
t-J model. In our calculations we use the same approximation, but previously obtain the effective vertices, which are 
perturbative. 

B. Types of pairing interaction. 

Generally, there are two different types of the hole-hole interaction in the t-J model. The first one is due to the 
spin-wave exchange, whereas the second is from minimization of the number of broken AF bonds by the placing of 
the holes at the nearest neighbor sites. We will consider them separately. 

The second type of interaction is usually introduced in the pure t-J model by adding projectors Pi = (1 — n 1 }) in 
J-term 

H J = jX)[(l-n?)SiS i (l-nJ)-Kni] > ( 36 ) 

id 

which project out the subspace of local states (spins), = h\hi is the operator of the hole number. It is evident that 
due to the h[{hi) and operator commutativity projection procedure is exact, i.e. there is no spin-spin interaction 
between the sites with the holes. Thus, the addition part to the t-J model can be written as 

&H t _j = J]T (-(r^ + n^SiSj + »f<S,S, - , (37) 

where summation run over bonds. Further treating this interaction term in the spin-wave approximation is straight- 
forward and yields 

AH t _j = -2J(1 - 2SX) ]T 7q /iL q 0k' +q 0k'/k + SHj , (38) 

k,k',q 

where term 5Hj, consisting of the two-magnon terms n h aa and n h n h aa, is presented in Appendix B. Hole attraction 
is enhanced by the zero-point fluctuations by the constant (— 2<5A ~ 0.16). Applying the transformation J25| ) to the 
Hamiltonian (|38| ) one can get the addition part of the effective Hamiltonian 

8H'° = J SViljL^^h ■ (39) 

k,k',q 

The expression of the second-order (in /k) addition to the bare vertex d38| ) is cumbersome and we present the general 
expression of SV^ 1 ^, in Appendix B. 
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Generally, there is an evident result for the n,nj interaction in the t — limit, namely, the groundstate of two 
holes is the bound state with the energy Ef, ~ — 0.58J. Moreover, there is the full degeneracy among the states of s 
[cos k x + cos k y ] , d [cos k x cosfcy], an d p [ sin k x , sin/c^] symmetries. The Ti n t ype of interaction ( |38| ) was intensively 
studied by the number of analytical [^2), and numerical techniques jio]], which established that increase of t 
(taking into account the spin-wave exchange interaction) leads to the rapid growth of Eb and disappearance of bound 
states. It was found that the largest critical value of t c is for the d symmetry of the bound state and that it is 
somewhere between t = (2 — 5) J, which is very close to the values of t proposed for the real CuC>2 planes. This fact 
has stimulated the developing of the idea of the superconductivity by the "preformed" pairs (bosons). 

Considered pairing interaction has nothing to do with the spin-fluctuation one, which was investigated in Refs. 
|f32f on the phenomenological basis and in Refs. J3(| using the RPA for the Hubbard model. An essential addition 
to the studying of the spin-wave exchange (Fig. pp) in the t-J and Hubbard models has been done in the work by 
Frenkel and Hanke pq| , where authors found that the exchange by the long-range (small momentum transfer) spin 
wave leads to the dipolar interaction between holes, which can be attractive or repulsive depending on the relative 
location of the holes. In the later work by Kuchiev and Sushkov |3(| this problem was independently studied in great 
details and several interesting features of the system were found. First of all, neglecting the retardation effect and 
finite size of the Brillouin zone for the two-hole problem one can obtain the Schrodinger equation with the potential 
~ (x 2 — y 2 )/r 4 , which leads to the "fall to the center" effect and to the infinite number of bound states. It was also 
found that the actual deepness of the bound states is very sensitive to the curvature of the hole band along the MBZ 
boundary and that for the equal masses in directions along and perpendicular to MBZ boundary there are no bound 
states at all. This effect is explained by the effective lowering of the dimensionality of the system when the mass along 
MBZ boundary tends to infinity, that is led to the strengthening of the pairing interaction. In Ref. |3(| only d and g 
[(cosk x — cos k y ) sink x sinfc y ] states were found to exist, and the lowest one was the d-state. This confirms a general 
statement | ]53[ that in AF state one-magnon exchange leads to the repulsion in s-wave and attraction in d-wave state. 

Note, that due to the long-range nature of the interaction and the large size of the proposed bound states it is not 
obvious how cluster analysis can directly reveal them. 



C. Bethe-Salpeter equation. 



Thus, one obtains the Hamiltonian (30) with the hole-magnon (|33| ) and "contact" hole- hole (|35|), ( p9|) interactions. 
Since we turned to an effective Hamiltonian by canonical transformation ( p5| ) the short-range spin-wave exchange 
(Fig. 0b) is now included in the "contact" interaction, which does not contain the retardation. Both interactions, 
according to above discussion, can lead to the pairing. As it was noted, the correct account for the retardation effect 
in the spin-wave exchange diagram is important, so let us consider this problem in details. 

The systematical method to search the bound states is to look for the two-particle Green function poles in the 
scattering channel as the function of the particles summary energy in the system of their center of inertia [pif . 
Corresponding integral equation for two holes with the total momentum P = is presented in Fig. ^| in a diagrammatic 
form. The standard way of solving of such a type of equation with the non-retarded "compact" bare vertex T° is 
given in Appendix C. 

In our case the " compact" vertex T° consists of two parts (see Fig. ||) and one has to include the magnon propagator 
into expression for the spin-wave exchange vertex. A natural assuming that the two-particle Green function has no 
singularities as the function of the difference of the energies of incoming particles provides somewhat different way of 
solving the Bethe-Salpeter problem. Details are also given in Appendix C. 

Resulting parametric equation on the bound state energy E of the Bethe-Salpeter type for the problem with two 
vertices (Fig. |^), which we will solve looks like 



^ E ) = E^2E^ 
p 



-2(fk ik+p M k , k+p ) 



Vj° „ + sv/ 9 



E - E p - £ k - UJq k,-k,q-T k,-k,q 

where q = k ± q for the exchange (direct) parts of V* 9 (|34|) and SV* 9 (|39|) . 



V(P,£), (40) 



D. Results. 

Thus, having in hands the vertices (^), ([35]), (|3^) and equation ([l0|) one can hope to obtain reliable results for the 
bound states in the t-J model. Moreover, since we have done the separation of the scales in the momentum space 
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and consider the " minimization of the number of broken bonds" interaction separately, one can demonstrate the role 
of each type of interaction in forming of the bound states. 

Shortly our results are as follows. The single bound state of d symmetry (d x 2_ y 2) is exist in all region of < 
(t/J) < 5. The states of other symmetries (s, p) at (t/J) > 0.2 were not found. The main thesis of this work is 
that the interplay of both interactions, which tend to <i-wave pairing, namely the short-range J-term (39) and the 
long-range spin- wave exchange (|3^), is important for the formation of this state. Specifically there are no bound states 
from the J-term alone for (t/J) > 2.1. The spin-wave exchange itself provide the shallow bound state, but the actual 
magnitude of the mass along the MBZ boundary almost pushes out this state in the continuous spectra. Nevertheless, 
putting these interactions together we get the bound state with the energy of the two orders of magnitude deeper 
than one found in Ref. |3(J. It was claimed in Ref. ||(| that due to the 1/r 2 nature of the attractive potential accurate 
account for the "cutoff" parameter (finite size of the MBZ) and the retardation effect is required. We have taken into 
account both effects and have done calculations for the lowest bound state. Deviation of the model from the pure t-J 
one, i.e. taking into account next hopping integrals, can lead to the different physics of the system. 

Let us discuss the other symmetry states. As it was noted above at the limit t = the J-term bound states of d, 
s, and p symmetry are not distinguished in energy, but already at (t/J) = 0.1 s- and p-states are two times higher 
than the d-state and they are disappeared at (t/J) = 0.2. 

Considering terms in equation in the bound state energy ( [lo|) separately and all together, we have obtained results 
for the ci-wave pairing state shown in Fig. [l0[ The dashed curve corresponds to taking into account the transformed 
J-term of interaction (|3^) alone (SV^ 9 _ k ). Obtained critical value of 2 for disappearing this, short-range in nature, 
state t c = 2.1 J is in excellent agreement with the variational approach |42]| , and in the good one with the finite-cluster 
calculations |28| , |Q , and other _approaches |l3| . Dashed-dotted curve corresponds to the long-range bound state due 
to the first two terms in Eq. (|f|). According to Ref. |36| this state should have small negative energy. We found 
that if one takes into account the second-order correction to the single-hole energy due to the rest of the hole-magnon 
vertex (J3G) it will lead to the slightly smaller mass along MBZ boundary and then to the very small binding energy 
~ — (10 — 10 — 4 )t for the long-range state. Actually, corrected value of the mass almost pushed up the bound state 
in the continuous spectra. As we noted in Sec. IV resulting masses coincide with the SCBA ones. The solid curve 
is our final result for the d-wave bound state energy in the t-J model. Bound state energy for (t/J) = 3 equals to 
AE = E — 2£k = — 0.0222;, which is of the two orders of magnitude deeper than it was obtained earlier |3(|. Thus we 
obtained the strong enhancement of the coupling effect because of the interplay of both types of pairing interaction. 

Note, that the "contact" part of the spin- wave exchange interaction (V? 9 ) plays the minor role in such a strong 
effect, namely ignoring it in Eq. ([Io|) yields the energy —0.012, which is only two times smaller than the result of the 
integral effect. 

It is useful to consider the structure of the wave functions of the two hole bound states in the k-space. Figure 
p"T| a shows the wave function for (t/J) = 1. It is simply "bare" short-range ipk ~ (cosfc^ — cosfcj,) with the small 
addition of the higher harmonics. Figures [Tl|b,c,d show wave functions for the (b) long-range state, (c) short-range 
one, and (d) resulting wave function, all for (t/J) = 2. Whereas the long-range bound state Ob is well localized 
near the band minima, the short-range one |Tl| c looks like something average between Figs, [Ha and [TT]b. It is not 
easily understandable, but one can say that in canonically transformed vertex 8V^ 9 ([59]) except the simple 7 q term 
there are some second-order terms with the same kinematic structure as in the V? 9 vertex, and since at (t/J) = 2 
this state is shallow the higher harmonics with this structure become important. Resulting wave function |ll| d reveals 
the features of the previous states. We believe that our results are very close to exact ones in all region of primary 
interest 1 < (t/J) < 5. Lowering of the bound state energy for the (t/J) = 5 manifests the worse accuracy of our 
approach for the very large t. 

The next problem, which can be touched, is the influence of the next-nearest hopping terms (2'-terms) on the bound 
states. Evidently small t' leads to the perturbative addition to the hole dispersion SE^ = 4t' cos(k x ) cos(k y ) and for 
the positive value of t' it makes the band more flat in (tt/2, 7r/2)— >(tt, 0) direction. According to the above discussion 
it also strongly enhances the long-range interaction and makes the d-wave bound state much deeper. For instance, for 
the absolutely flat band (m\\ = oo) at (t/J) = 3 energy of the bound state is E = —0.1652. Note, that the neglecting 
of the short-range interaction (J-term) provides result of the next order of smallness E = —0.0232. Further increase 
of 2' will lead to the shift of the minima to the (0, ±7r), (±7r, 0) points and will make the bound state shallower. For 
some region of 2' > the long-range bound state of g symmetry becomes possible. Its wave function obeys the same 
symmetry as [(cosfc x — cos k y ) sin k x sin/cj,], i.e changes the sign in MBZ eight times. Since there is no short-range 
interaction for this state, the energy associated with it is very small. 

It is well established by now, that for the real Cu02 compounds 2' has the negative sign and the total effect from 
2'-terms is the fully isotropic dispersion near the band minima |5^ , |56]| . Note, that it is the main effect from such 
terms and one can neglect their contribution to the effective vertices up to the rather large values of 2'. Joining this 
statement together with the claimed sensitivity of the bound state to the anisotropy of the band one can suggest that 
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there are no bound states in the realistic model of the CuC>2 plane. We have studied the problem of the critical value 
of t' and found that the difference A = E(n,o) ~ E(n/2,n/2) i which is directly connected with the inverse mass along 
MBZ boundary, at which the bound state energy tends to zero is equal to 0.22i for (t/J) — 3. Comparing it with 
the bare t-J model value of A = 0.13t one can get the critical t' c ~ 0.3J, which is much lower than the realistic value 

t'eff ~ 1-5J. 

Turning back to the simple t-J model, one can say that, at the first glance, obtained result strongly supports the 
idea of the t-J superconductivity by the condensate of the "preformed" pairs. In our opinion such a picture is not so 
obvious, since the existence of the long-range spin wave excitations is essential for the bound state formation, whereas 
the long-range order is unstable under very small |QJ5| finite concentration of holes. Therefore, for the clearing this 
subject one has to solve the pairing and stable spin state problem selfconsistently. 



VI. CONCLUSION 



We conclude by summarizing our results. 

We have put forward the canonical transformation of the t-J model Hamiltonian using an analogy with the lattice 
polaron problem as well as some evident ideas based on the known single-hole properties in the AF background. We 
have shown that the rather wide type of the transformations, which has the some type of the 1/z expansion in the 
basis, allow to extend the region of the analytical treatment of the problems up to t/J ~ 5 with an appropriate 
accuracy. Generally, powerful method applied provided us the straight way to the formulating of the quasiparticle 
Hamiltonian, which includes the free energy terms and all essential interactions. 

The proof of the correctness of the simple form of the effective quasiparticle Hamiltonian has been done. Results 
for the single-hole bottom of the band, bandwidth, and other properties have been compared with ones of the SCBA 
calculations and a remarkable agreement has been found. It is supported the idea that the "canonically transformed" 
quasiparticles have the properties, which are close to ones of exact t-J model quasiparticles. Although the bandwidth 
and components of magnetic polaron wave function have discrepancies with the exact ones, which are grown for the 
larger (t/J) > 5, the bottom of the hole band remains very close to an exact one, confirming the basic idea that the 
effect of the polaron well formation (of the order of t) is taken into account by our transformation. 

Using the obtained Hamiltonian we have performed the study of the two-hole problem. The hole-hole interactions 
of the different nature have been considered separately and all together. Rather deep bound state of d-wave symmetry 
originated from the interplay of two types of the pairing interactions has been found. Retardation effect for the long- 
range spin- wave exchange has been carefully taken into account. Other possible symmetries of the bound state wave 
function have been studied as well. The main effect of the so called t'-terms has been investigated and the critical 
value of t' c for the bound state existence has been found. 

Since we have used the presence of the AF long-range order as a foundation of the setting up the problem and the 
long-range interaction was found to be essential for the bound state formation, so the direct relation of the considered 
two- hole problem to the case of finite hole doping of the real Cu02 plane is unclear. We have briefly discussed the 
questions, which remain to be resolved. 
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APPENDIX A: 



The series of commutators of the type ( p5| ) with generator of the transformation (^) provide the general expression 
for the hole energy and hole-magnon formfactor (ft) 



E k = 2tJ2 (^k,qMk,q 



gfWk.qA'k.q 



5T<5Wk,qMk,q 



2j r 

t 



1 n 2 , , 
STPk^q^' 



q I" J[Wk,qAtk,q + gf<Wk,qMk,q 



(Al) 



13 



2J r 



F k ,d = 1+1 5T W k ,q + JT<5W k ,q + — Mk,q^q + |f W k ,q + ^<W k ,q /M k 



with 



W k ,q = 3/i k , q Vfc-q - Vk - M k>q C/ k _ q + f/ k , VK k ,q = ^k,q 



3 Vk-n-F k -U>q f/k-q + t/. 



Wk,q - Mk.q 15(f/ k Vk - C/k-q^k-q) + 5(C/ k _qF k - f/ k Vk_q) + M k ,q + + 6i7 k _ q 17 k , 



A i k,q|5([/k-q'A^k,q'A i k,q' — £4:-q-q' Aik-q.q< /lk-q,q' ) + 8(V k _ q _ q '^^_ q q , — Vfc_q' A^q' ) 



Mc.qi ^k-q'Mk.q' + ^k-q-q'^k-c 



(A2) 



<*W k , q = /i k ,q( 15(t/ k t/£ - C/ k -qC/ k J -q) + 5(tf k _ q £^ ~ ^k t^_ q ) + W q ( [7*_ q + ^ + 6 (7 k _ q[ / k ) 



£q< 



5(t/k-q'A i k,q' ti 'q' _ ^k-q-q' Mk-q.q'^q' ) + K^ik-q-q'Mlc-q.q' _ ^k-q'Mk,q') 
+W q {?7k-q'Mk,q' + ^k-q-q'Aifc-q.q' } 



^ k T q 



where 



Vk = Y^ M k,q'Mk,q', t/ k = ^k,q'' ^ k = ^ /4,q' W q' , 
q' q' q' 



(A3) 



and 8W£q, SW^^ are the terms from the Trugman processes. They were found to play the negligible role for the 
energy and vertex. 



Explicit expression of the energy up to the fourth order over the TP for the generator in the form (25) is 

J 



£k * 2i (A4 + | [3/^ - 4/^72 - /k^if] + - [f& + i (3/^k 3 ~ 4/ k 4 / k /k - 

With F k ' = £ (A/2 , q /k-q/ k -q) , F£=J2 (^q^-q/k-q) , 



(A4) 



*2 



E ( M k,q^k-q/L q ) , ( M k >q ^-q/ k 2 _ q Wq) , 



Y± = E ( M k + q,q) > ?k = E (M2 +q , q 6 



where J k and J k are defined by Eq. ( |27| ) 



Hole-hole vertex V^, q = V c h x h (k, kVq) + V^(k, k', q) is 

V^(k,k',q) = ( A/ k , q Aik<+q,q + A/k-q,-qMk',-q 



1 + ^T^k.k'.q 



(A5) 



2.7 



Wq I /Uk,q/^k'+q.q + /^k-q,-qMk',-q 



2! + 4! ^k.k'.q 



4 I A^qMk'+q.q — Mk-q.-q^k'.-q 



9 J 



with 



P k , k ',q = -3(C/ k _ q + U k ) - (tf k ,_ q + ) , 



and 



KttrO 4 ' ^ ! C l) — ^k, k ',q + ^lc-q.k'+q^-q + ^k,k\q + ^lc-q,k'+q.-c 



(A6) 



^Ic.lc'.q — E Mk'+q,q-q' 



A* k '+q',q'(^k-q%q-q'/^ k ,q' _ M k) q' Mk-q' ,q-q' ) 
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Ti 



2 

k,k',q 



22 Mk'+q.q'Mk-q+q'.q' 

q' 



Mk',— q' (-^k— q',q— q'/^k— q',— q' ~t~ -^k— q',— q'Mk— q',q— q' 

^k,q-q'Mk'+q-q',q-q' + J"k-q+q',-q+q'A t k',-q+q' 

+A/k'+q— q'.q— q'Mk.q— q' + ,— q+q'/-*k— q+q',— q+q' 

^""(k, q, q') = M k _ q , q /^ k , q - M k , q Mk- q ,q' + — (w q - w q /)^i k _ q!q /^ k , q + higher order terms 

V^ aa (k, q, q') = Mlc^q/ifc-q-q/^q' + ^qA/lc-q-q'^q/ — M k _ q ' j q / U k _ q ' ! _ q / — ^k-q',q^k-q',-q 

+ j(w q + w q /)[/ik,qMk-q-q',-q' - Mk-q', q Mk-q',-q'] + higher order terms 



Two-magnon vertices are 



(A7) 



APPENDIX B: 



Hole-two-magnon addition SHj (38) originated from the projectors applied to J-term has the form 

SHj = J { V qL> (4-q-q-q^L+Q^/ika^ + H.c.) + 2^ Q 7 ^ q ,/ lk _ Q - q+q ^L+Q^/ l ka q ^ } (Bl) 

- J E {^q J q^4-q-q^ka q a q ,+H.C.)+2y4^tq + q^ka qaq j , 



k,q,q' 

where vertices are given by 



7Q(""qWq' + WqUq') + \ [iQ+q + 7Q+q']( M q U q' + "q^q') > 
^q.q' = 7Q("qUq' + WqWq') + \ [7Q+q + YQ-q']0qV + «q"q') » 



Q,q,q' 

J, 2 



q,q' 

J,2 



[1 + 7q+q']("q"q' + VqMq') + frq + 7q'](w q Wq' + Wqfq') 



Kj.'q' = I 1 + 7q-q'](w q Uq< + l> q U q ,) + [ 7q + 7 q '](u q V q , + W q U q /) , 



(B2) 



where u q , w q are the parameters of Bogolubov transformation, C/ k defined in (A3). 

Applying the transformation (25) to the Hamiltonian (Bl) one can get the additional part to the effective interaction 

© ^k.k'.q = 5V dir^ k '> q) + 5V ei 9 (K k ', q) with the direct part 

«^(k,k',q) = - 7q (l - 2SX)[l - C/ k - (7 k _ q ] (B3) 

+4^Mk,q'/k-q' |<5A7 k _q/(7 q u q ' + 5 7q -q'Uq') + 5 u q'-Pk-q' : q} + [k, q => k q, -q] 
q' 

+4^/Xk',q'/k'-q' ||5AUq'7k'_ q ' 7 q +q ' + iu q >.Fk'-q',-q + Uq'^k'-q'.-q-q' | + [k', q => k' + q, -q] 

q' 

" Mk'+q.q-q' {/^k'+q'.q' V^q-q/.q' + Mk'.-q' Kj-q' ,-q' } + [ k '> q ^ k' + q, -q] , 

q' 

and exchange part 

SV-//( k ,k',q) = 32/ k / u , {i<5A 2 7k7 k'7q + <5A 7k F k ,,_ q + ±D llk ,F k<(l + i£> 2 F k , k , +q + i£» 3 F k> _ k , +q } (B4) 

-2 {/Xk'+q,q$k,q + $ k , q } + (1 - 2<5A) 7q+q'Mk',q'Mk+q',q' 

q' 

+[k,k',q=^k-q,k / +q,-q] , 
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where functions 

F kl ,k 2 = ^llkiTka + i-D27ki-k 2 + i-D37ki+k 2 , (B5) 
*k, q + *k, q = 4/k {^A 7k (7 q w q + « q ) + (i + <SA)u q7k + « q £ k , q } + [k, q => k q, -q] , 

with the numbers 

SX = ^(^ q + 7q«qUq) = —0.0790 , D\ = ^(7qt>q"q + u q cos q x cos q y ) = —0.1301 , (B6) 

q q 

D 2 = ^2 u q(! - cos Qx cos q y ) = 0.2628 , D 2 = ^ u q (2 cos 2 tfe - 1 - cos q x cos q y ) = -0.0077 . (B7) 
q q 



APPENDIX C: 



For two holes with the total momentum P = one can write the following integral equation 

f(k f , -k g , k' f ,-k' g ) = r°(k f -k g , k'f, —kg) + r °( fc /' -k 3 ,Pf,-P g )G( Pf )G(- P g)T(p f , -pg, k' f ~k' g ) (Cl) 



pf 



where we introduced four-momentum notations kf = (k, e/), — k g = (— k, e g ), pf — (p, e'^), —p g = (— k, e' g ) with 
momenta k, p and frequencies e "f( g )- G(p) = l/(e — £7 k + i£) is the single-hole Green function. This equation 



is equivalent to the graphical equality shown in Fig. Near the pole r° <C T and hence the first term in Eq. (Cl) 
can be neglected. Then, one can see that T dependence on outgoing four-momenta k, kl is the parametric one, i.e. it 
is not determined by equation, and they can be omitted. Let us also introduce E — ey + e g , Ae = (e/ — e 9 )/2, and 
," _ , " 



Ae" = (e" - e")/2. Thus, 



f (k, £7, Ae) = Y, r °( k ' P. E > Ae ' Ae ")G( P , E/2 + Ae")G(-p, £7/2 - Ae")f (p, £7, Ae") . (C2) 

p,Ac" 

When r° has no frequency dependence ("static" interaction), 

r°(k,p,£,Ae,Ae") =t/(k,p) , (C3) 
it is natural to change the variable GGT = \ and get 

X (k, £7, Ae) = G(k, E/2 + Ae)G(-k, £7/2 - Ae) £ U(k, p) f d(Ae")x(p, E, Ae") . (C4) 
Additional integration of both sides over Ae gives the Schrodinger equation: 

V>(k, £7) = - Y U ( k > pMp, , (C5) 

p 

with ^(k, £7) = J (i(Ae)x(k, £7, Ae), which has the sense of the bound-state wave function. 

In our case the " compact" vertex T° consists of two parts (see Fig. ||) and one has to include the magnon propagator 
into expression for the spin- wave exchange vertex 



r?(k,p,Ae,Ae'' 



^,P^k,-p * x pl k 



e - e" — Wk+p + iS e" — e — cj_ k _ p + iS 



(C6) 



where Vk, P is the vertex of the type Eq. (|33|), e — e" = Ae — Ae", k + p = q. Both negative sign and q = k + p are 
due to exchange character of the diagram (Fig. ||). Thus, there are three Ae"-dependent denominators in the integral 
Eq. dC2] ) and the above change of variables is impossible. 

It is natural to assume that since one is looking for the poles of the two-particle Green function as the function of £7, 
r has no singularities as the function of the difference of the energies of incoming particles. Therefore, the integration 
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over Ae" in Eq. (|C2[) can be done and it is determined by the poles of G(p,E/2 + Ae"), G(p,E/2 — Ae"), and 
r?(k,p, Ae, Ae") (jcff. Their poles are: Ae" = (E p -E/2)-i6, Ae" = -(E p - E/2U-iS, and Ae" = ±(Ae- w q ) ±iS, 
respectively. +(— ) in the last pole corresponds to the first (second) term in Eq. dCq). The integration gives 



f(k,£,Ae) = £ 



E-2E n 



T(p,E, (E p -E/2)) 



T(p,E,-(E p -E/2)) 



Ae - {E p - E/2) -LU^ + iS -Ae - (E p - E/2) -tu^ + iS 



(C7) 



The further way is close to the usual one. Multiplying both sides of Eq. (C7) by the external incoming Green functions 
one can integrate over Ae, using the evident parity of T on Ae 



E - 2E U E - 2E k ^ E - E p - E k - w q X E- 2E p 



(C8) 



Changing T(k, E)/(E - 2£' k ) = t/>(k, E) one obtains 



^(k,E) = 



-2Vic, P VlV-p 
E- 1 — ^ — ' E — Ep — E\t — cjq 



<Kp,£) 



(C9) 



Evidently, the "usual" Bethe-Salpeter equation (C5) can be obtained by the same way. Surprisingly, this result (|C9| ) 
coincides exactly with one obtained in Ref. [Bm using Rayleigh-Schrodinger perturbation theory. 
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Figure captions 

FIG. 1. Single-hole energy for the Ising limit. Bold solid curve is an exact result in the spin-wave approximation. Dashed 
curve is an exact result for the large t ^> J limit. Solid curves (1) and (2) are the canonical transformation results up to the 
fourth and sixth order of transformation, respectively. 

FIG. 2. Weights of the components of the magnetic polaron wave function (Ising limit). Bold curves are for an exact 
solution, light curves are for the canonical transformation results (fifth order). Solid curves correspond to the weights of 
the bare hole, dashed-double-dotted curves - hole+1 magnon, dashed curve - hole+2 magnons, long-dashed curve - hole+3 
magnons. 

FIG. 3. Bottom of the hole band. Solid curve is our result (the sixth order of the transformation), dashed curve is the 
SCBA result. 

FIG. 4. Width of the hole band. Solid curve is our result (the sixth order of the transformation), dashed curve is the SCBA 
result. 
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FIG. 5. Weights of the components of the magnetic polaron wave function (Neel case). Bold curves are for an exact solution 
(SCBA), light curves are for the canonical transformation results (the fourth order). Solid curves correspond to the weights of 
the bare hole, dashed curves - hole+1 magnon, dashed-dotted curve - hole+2 magnons. 

FIG. 6. t/J dependence of the formfactor -Fk, q . Solid curve -k = (7r/2,7r/2), q = (0,0), dashed curve - k = (7r/2,7r/2), 
q = (tt/2, tt/2), dashed-double-dotted curve - k = (0, 0), q = (0, 0). 

FIG. 7. Schematic view of the scattering diagrams: (a) // — > //, (b) fg — > gf of exchange type, (c) fg — * fg of direct type, 
(d) the first diagram of the fg — > fg scattering beyond the ladder approximation. Here, the wavy lines denote the interaction, 
which originates from the magnon exchange (t-term), and the point in the diagram (c) denotes the vertex, which is from the 
nearest neighbor interaction (J-term). 

FIG. 8. Graphical identity for an exact vertex T(kf, —k g , k't, —k' g ) for the fg scattering in the ladder approximation. Black 
circle denotes F(fc/, — k g , k't, —k' g ), empty circle denotes a "compact" vertex r°(fc/, — k g , k't, —k'g). 

FIG. 9. Structure of the "compact" vertex F°(kf, —k g , k't, —k' g ) (empty circle). Here, the wavy line denotes the interaction 
due to the long-range spin- wave exchange (p3|), and the point denotes all diagrams, which do not contain the retardation (^B[), 

©■ 

FIG. 10. Results for the energy of the d-wave pairing state. Dashed curve corresponds the short-range bound state, 
dashed-dotted curve corresponds to the long-range one, and solid curve corresponds to the combined effect of both types of the 
state. 

FIG. 11. Wave functions of the two-hole bound states: (a) (t/J) — 1, (b) (t/J) = 2, long-range state (c) (t/J) = 2, 
short-range one, (d) (t/J) = 2, wave function of the mixed state. 
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